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Abstract 

A study is presented for the non linear evolution of a self gravitating 
distribution of matter coupled to a massless scalar field. The characteris- 
tic formulation for numerical relativity is used to follow the evolution by 
a sequence of light cones open to the future. Bondian frames are used to 
endow physical meaning to the matter variables and to the massless scalar 
field. Asymptotic approaches to the origin and to infinity are achieved; 
at the boundary surface interior and exterior solutions are matched guar- 
anteeing the Darmois-Lichnerowicz conditions. To show how the scheme 
works some numerical models are discussed. We exemplify evolving scalar 
waves on the following fixed backgrounds: A) an atmosphere between 
the boundary surface of an incompressible mixtured fluid and inflnity; 
B) a polytropic distribution matched to a Schwarzschild exterior; C) a 
Schwarzschild-Schwarzschild spacetime. The conservation of energy, the 
Newman-Penrose constant preservation and other expected features are 
observed. 

Key words: Characteristic Formulation; Matter Evolution; Einstein- 
Klein-Gordon System 

1 Introduction 

Numerical relativity has reached high sophistication levels to advance in the 
study of realistic solutions to the Einstein equations. [1 . Particularly, to sim- 
ulate gravitational radiation from collapsing sources ^2] and binary black hole 
mergers [3] . Almost all the investigations have been done in the ADM 3 + 1 for- 
mulation [4j , although the characteristic formulation offers a myriad of valuable 
advantages 0. 

This work is motivated by the possibility of simulating gravitational radia- 
tion from axisymmetric matter sources, as part of a longer project [Sj, [Zj- We 
met that purpose beginning with the most simplified and concomitant problem 
under spherical symmetry, that is, the self-gravitating massless scalar field. In 
the same vein we deal with the problem of matter coupled to radiation which 
eventually will lead us to get gravitational signals from bounded sources. As 
far as we can see, this problem stands alone as an important one in the field. 
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Here we report a study in the aforementioned direction, limiting ourselves to 
spherical symmetry to follow the evolution of a massless scalar field interacting 
with a perfect fluid distribution of matter. This model problem offers a number 
of advantageous computational and geometrical features. It is well known that 
the scalar field mimics gravitational radiation and it has been used to study the 
global properties of the spacetime, black hole threshold and radiative signals 

i, El, [10], [a, [la. 

The framework for computing a complete spacetime within the characteristic 
approach has been laid out by Tamburino and Winicour [13] and explicitly in 
|14j , where is explored the production of gravitational waves for axially symmet- 
ric distribution of matter (for a complete review see [5] and references therein). 
Some years later, fundamental studies in ID and 3D [5], [TS], [H] were followed 
by a 3D code able to treat matter by Bishop et al. ^Tj}.. A relatively simpler 
system is a spherically symmetric distribution of a perfect fluid coupled with 
scalar radiation [TS]. For other symmetries other than spherical, the mathemat- 
ical problem has a very similar structure [2] . In these contexts we can take full 
advantage of the characteristic approach to treat matter and radiation. The 
first example of the use of characteristic numerical relativity for the study of 
dynamical neutron star spacetime, collapse and radiative signals was reported 
in [TB] . Few investigations consider the scalar field interaction with fluid stellar 
distributions [IH], HO]- 

We have recently discovered an unexpected unity in the treatment of mat- 
ter in numerical relativity [2T], using explicit Bondian observers |22J- These 
observers offer an Eulerian (noncomoving) description (global) with the spirit 
of Lagrangian observers (local and comoving). Following this line also we re- 
ported a disclosure of a central equation of state (CEoS), which is unique for 
all evolutions; it emerges as a conserved quantity from the field equations [53]. 

Authors commonly refer to the scalar field as a model of matter distribution 
that simplifies the treatment of the hydrodynamic issues. This approach has 
been useful to study non linear physics and asymptotic behaviors, especially 
for central regions [23]. Scalar field models have been extended to realistic 
situations such as gravitational radiation [l5]. However, cases where the scalar 
field is coupled with radiation require a different approach due to both the 
confinement effect of matter and the dispersive nature of radiation. When a 
scalar field is coupled to matter it can be easily interpreted as an anisotropic 
fluid if we use Bondian frames explicitly. 

We develop a numerical framework to deal with matter coupled to scalar 
radiation. We perform a detailed study of the central world line at r = 0, 
leading to a conformally flat spacetime in that region. We assume that the radial 
dependence of the geometrical and physical variables keep the same dependence 
as the static variables near the coordinate-origin. 

We formulate in this work the characteristic evolution in terms of Bondian 
observers [3TJ [53] . The Eulerian formulations of numerical relativity [TS] , [55] 
actually use Bondian observers in the mathematical treatment of matter. We 
also match the interior solution with the exterior one in a clear and precise 
treatment of the boundary distribution of matter without the use of artificial 
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atmosphere. Infinity is treated as usual in literature. 

In what follows we write the field equations for Bondian observers when mat- 
ter is coupled to scalar radiation, which makes the fluid manifestly anisotropic. 
Section 3 is devoted to the regularization and matching. In section 4 we show 
how the scheme works by means of numerical test models. Finally we summarize 
with some remarks in section 5. 



2 Field equations for Bondian frames 

Bondi's metric in the spherical form reads [22] 

ds^ = e^P(-du^ + 2dudr) - r^{de^ + sin 9^ d(t)^), (1) 



where /3 = I3{u, r) and V ~ V{u, r). In these coordinates the components of the 
energy-momentum tensor are distinguished by a bar. In spherical symmetry 
there exists a well defined notion of quasilocal energy, the Misner-Sharp mass 
function, m{u, r) [55] introduced by means of 

m= i(r-ye-2'5), (2) 

which measures the energy content in the sphere of radius r and it reduces to 
the Arnowitt-Deser-Misner and Bondi masses in the appropriate limits. 

Consider a stress-energy tensor for a perfect fluid and a massless scalar field 

%,^f^l+fl. (3) 

One can follow the Tamburino-Winicour formalism [T3], in particular as ap- 
plied in regular spacetimes, where the foliation of light cones emanates from a 
freely falling central observer [I1]-[IH]. But following Bondi, local Minkowski 
coordinates {t,x,y,z) are introduced by 

dt = e2'5(l - 2m/r)i/2^u-f (1 - 2m/r)-i/2dr, (4a) 

dx = (1 - 2rn/r)-^^^dr, (4b) 

dy = rdO, (4c) 

dz — r sin 9dip. (4d) 

Denoting the Minkowski components of the energy-momentum tensor by a caret 
we have 

Too = fooe^^{l-2m/r), (5a) 
Toi = (foo+foi)e2^, (5^^ 
Tn = (l-2m/r)-i(foo + fn + 2foi), (5c) 

f 3 



T| = Tl^fi=fl (5d) 
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Next one assumes that for an observer moving relative to these coordinates 
with velocity w in the radial direction, the space contains an isotropic fluid with 
pressure p and energy density p. 

For this Bondian observer, the covariant energy-momentum tensor of matter 

is: 

/ p \ 
p 
p 
\ p J 

Then a Lorentz transformation readily shows that 



rpM 



(6) 



Too 
Til 



l+UJ 



p + puj-' 

1-^2 



(1 - 2m/r)-\p + p) 



l-OJ 
l+OJ 



-P- 



(7a) 

(7b) 

(7c) 

(7d) 



The energy-momentum tensor for the massless scalar field minimally coupled 
with gravity 

1 



(8) 



can be read by an observer at rest in the frame of Q. 

It can be shown that the Einstein-Klein-Gordon equations can be written 

as 



1 — oj^ 47rr 



r — 2rh 



(9) 



p-up ^ m,r 

' p - 



1 +0; 



47rr2 ' 



i-^(p + p) + +p* = (1 - 2m/r)^, 
1 + u! Zirr 



(10) 
(11) 



P + Pf = -^/3,„,e-2^' + -L(i„2m/r)(2/3.„ + 4/32^-/3.Jr) 

47r OTT 

+ 7r'-[3/3,r(l - 2m.^) - m.„] (12) 
oTrr 



and 



(13) 



where the comma denotes partial differentiation respect to the indicated coor- 
dinate, and the scalar energy flux e*, the scalar energy density p*, the scalar 
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radial pressure p*, the scalar tangential pressure pf , respectively are defined by 
e* = e-2^[e-2'3(l - 2m/r)-^<^\ - (14) 



9* = = (1 - 2m/r)$y 2, 



(15) 



(16) 



From this point of view the scalar field can be interpreted as a radiating and 
anisotropic fluid ^\ whose energy-momentum tensor can be written as 



(17) 



with u^u^ = 1, Z^/p = 0, x^Xfi = ^li if we identify the four velocity for an 
observer at rest in the frame of ([T]), the null and the space-like vectors as 



(18) 



(19) 



X^ = (l-2m/r)-V25i. 



(20) 



Note that Bondian observers can be purely Lagrangian when we deal only with 
radiation (21j . 

The conservation equation T|!^ = 0, or equations ([To|-(|T2|), lead us to the 
generalized Tolman-Oppenheimer-Volkoff (TOV) equation for aon static radia- 
tive situations 



P,r - e 
2 



-20 ( P + P \ . ( P+P 

\l-2rfi/rJ,u Vl-2m/r 



47rr(p + pf) + mjr^ 



r 



ip-p), 



where 



P = 



p_--ujp 

~ p-t^p 

1 + w 

are the named effective variables [55], 

From the field equation ^ is straightforward that 

^ = -Attt^Hp - p*)^ + (1 - 2m/r)e2^e*] 
au du 



where 



— =e2/5(i_2m/r)- , 

du 1 — cj 



(21) 

(22) 
(23) 

(24) 
(25) 
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is the matter velocity. Integrating the wave equation (13) we obtain 



Jo 



(26) 



which combined with ( 24 ) lead us to 



du du 



+ 27rr^(l - 2m/r) ( ) e^^p*, (27) 

1 — UJ J 



where 



1 r 

N = -e"2^ J e^^{l ~ 2m/r)$,,.dr. 



In absence of matter, p — p — 0, for the exterior region the field equations 
reduce to 

R^, = -8^$,^$,,, (28) 
or explicitly to the hypersurface equations 

P^r - 27rr$2^, (29) 
m = 27rr(r - 2m) $2^, (30) 



and to the wave equation (13) [8]. 



3 Regular izat ion and matching 

Some previous investigations consider regularization near r = [H], [^S], f^ . 
[30] . |31) . The conditions for the scalar field, as a matter model, do not necessar- 
ily apply to distributions of matter. Depending on gauge conditions each proce- 
dure to get regular spacetimes may be cumbersome and tricky, even in vacuum. 
Initially regular spacetimes can eventually develop singularities [21], [22]. We 
show a simple way to construct regular spacetimes, near the coordinate-origin, 
when the inner spacetime corresponds to a spherical distribution of baryonic 
matter coupled to a massless scalar field. To construct regular and general 
enough spacetimes, which eventually recover equilibria, collapse, form singular- 
ities and horizons, we do an asymptotic study close to the special regions: r = 
and r — >■ oo. The treatment is basically the same for these two zones, that is, 
power expansions of r and r^^, respectively. For r = R{u), the boundary sur- 
face, the Darmois-Lichnerowicz [33 j .|34j conditions are guaranteed to match the 
interior and exterior solutions on a moving boundary. Integrating from r = 
no additional conditions are required at the surface to describe its evolution. 
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3.1 Close to the origin 

Consider the following asymptotic expansions for the metric functions near r 
0, which represents a regularly and conformally flat spacetime [53]: 



TO = TO3(u)r^ + 0{r^), 

and for the escalar field 



(31) 
(32) 

(33) 



Using the field equations we get the physical variables expansion as a function 
of r as showed in Table I. 



o 


np 


np 


w 





|to3 - ^n^'i 







1 


Wi(;92-7I"$^)+7r<I>i-I>2 


wi(/32 - n^l) - 27r$i$2 





Table 1: Coefficients for the physical variables expansion as a power of r near 
the center. 



From the evolution equation ( 13 ) we get 

2/30^^0 _ 



du 



$1, 



and 



3 

du 2 



From the field equation ( 12 ) 



du 



[7r$i$2 - 2wi(/32 - 7r$i)], 



and from ([9]) 



(34) 
(35) 

(36) 
(37) 



This last equation, together with expansions showed in Table I, is readily inte- 
grated to give 

Pc + 3pc = constant, (38) 

where Pc Po + pt and Pc^Po+pt- 

The lapse of the coordinate time du is related to the corresponding lapse of 
time dr measured by a central observer as 



dr = e^^°du. 



(39) 
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We have preference for the central time r to describe the studied system. The 
reason is mainly numerical: proceeding with the radial integration from r = 0, 
where the world line is therefore geodesic, in consistence with the conformally 
flat result as a consequence of regularity. Therefore, the replacements P — >■ /3— ,80 
and u ^ T left invariant the field equations (|9])-(13). The same situation will 



be analog for the matching surface, which behaves asymptotically as Vaidya or 
Schwarzschild, and for the asymptotically flat infinity. 

3.2 Matching at the surface 

Boundary conditions at the surface r = R{u) are needed for j3 and rh in order 
to perform radial integrations. We match the interior solution with the exterior 
at r = i? by means of the Darmois-Lichnerowicz conditions. These conditions 
are equivalent to the continuity of the functions /? and rfi across the boundary, 
and to the continuity of the spin coefficient [34] 

7=(l-2m/0/3,r-^-/3,„e-2'3. (40) 
Considering the expansion of (3 around the surface, we have 

where the superscript ± indicates the evaluation of the function aX r — R + Q 
or r = i? — 0, and B = /3{u, r — R). Therefore, 7+ = 7" leads us to 

u:j,^l-2RF^, (42) 

where F = 1 — 2M/i?, M — rhj^, the subscript R indicates that the quantity 
is evaluated at the surface, and (5^* — — represents the jump of the 
indicated function across the boundary. Observe that SR = 0, where over dot 
indicates derivative respect to time. Thus, we get 

Pr + ujuPb = {l + i^B)6p'^, (43) 

which leads us directly to 

PH=(\^)sp-. (44) 

If the scalar field gradient is continuous across the boundary r = R we have a 
pressureless surface, which can be expressed as pR = —ujrPr or equivalently 

m--2RFI3-+inR^p% 
= — — , n2 • (45) 

m^r — iT^R pJi 

When PR — 0, we solve as usual the indetermination taking the limit 

w_R = — lim ^ = — ^ , (46) 

■r-s._R p p r R 
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where 



-{p- 
r 

27rr2 



{Am - 2m^r - I) - p^, 



and 



P.r 



(47) 



(48) 



Two comments are in order here. First, ( |46[ ) can be obtained from the continuity 
of the spin coefficient gradient 7^ = 7". Second, exactly the same result is 



obtained from the field equations, that is, equation (461 proceed from (21 1 
evaluated at the surface ^28]. 

Once satisfied the matching conditions across the boundary r — R we need 
surface equations to follow the radius evolution and the exchange of energy on 



it. Evaluating (27) at the surface and defining 



and 



we get 



7^* = 4^i?V|, 



-2B 



dM 

du 



-AttNI 



(49) 



which clearly establishes the transfer of energy at the boundary of the matter 



distribution. Additionally, evaluating ( 25 ) at the surface we get 

,dR 



-2B ' 



du 



(50) 



which conforms together with (45) and (49) the system of equations at the 



surface for pji ^ 0. These equations determine completely the evolution at the 
surface distribution. 



3.3 Close to infinity 

This section is standard in literature but we include here a resume for the sake 
of completeness. Assuming that the scalar field has an asymptotic expansion 

i 



r 



the metric functions read 



2ttQI , 27rQi(7WQi-2Q2) 
r 



m = M{u) - + -^^^'-^ + o{r-^). 



(51) 

(52) 
(53) 
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The coefficient Qi is the scalar monopole moment and Q2 the Newman-Penrose 
constant. The asymptotic Bondi mass A4 also can be expressed globally as 



Ai{u) = 4:71 / r^pdr + 4-K I p'^ dr — m\j+ , 
Jo Jo 

and the scalar news function [5] 

1 r°° 

M{u) = / e2^(l - 2m/r)<S>,rdr. 

2 Jo 



The Bondi mass loss equation is 



e —r- = -47r7V^^. 
au 



(54) 



(55) 



(56) 



With these definitions, the total radiated energy can be established 

A = M{u)-M{uo) = -J{u), (57) 



where 



(58) 



It is easy to check from ( 54 ) that the most general Killing propagator to get the 
energy conservation in the spherical context is 

= e~^'^S^. (59) 



As a matter of fact, we get (59) when the Linkage [T^], a finite representation 
of the Bondi-Metzner-Sachs asymptotic group. 



(60) 



is compared with (54). 

Up to now the system has been described with enough generality to proceed 
with the hydrodynamical solver developing, which was not planted as a goal in 
this paper. However to show how the scheme works some numerical models are 
discussed. 



4 Models 

4.1 Ghost scalar fluid 

This model exemplifies the junction conditions, the scalar field energetics and 
the preservation of the Newman-Penrose constant. 

Consider a mixture of two components fiuid of pressure p + p* — > P and 
energy density p + — > p which is incompressible and remains static if a layout 
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outer scalar field guarantee the Darmois-Lichnerowiczs conditions. Clearly from 
(45) these conditions reduce to 



m- - 2RFI3^^ + 47ri?VI = 0, (61) 

which is satisfied if 

FI3+ = 2ttRpr. (62) 

Under this scenario uj = everywhere, and derivatives respect to the timelike 

coordinate dropping to zero. The set of equations 12 1 simplify considerably 
to get the well known interior solution [55] 

m = M{r/Rf, (63) 



=2/3 




where 
and 



P 



pl/2 _ ^1/2 
^1/2 _ 3^1/2 

[l-il~F)ir/Rf 
3(1 - F) 3Af 



(64) 
(65) 
(66) 
(67) 



Observe that (3r — 0, F(3'^ = AirRpn = 3M/ R^ and pn. = 0. In this way we have 
constructed a "fixed" background of a incompressible perfect fluid coupled to an 
exterior scalar field which remains "frozen" in some hydrodynamic characteristic 
time scale. 



$ 


Family 


Parameters 


X/{R + r) 


(a) 


A, R 


X{ra-rnrb-r)y[in-r,)/2r 


(b) 




Aexp (r — tqY ju'^ 


(c) 


A, (7, ro 



Table 2: Exterior/Initial data set for the massless scalar field 



Three different families of exterior /initial scalar fields are showed in Table 
II. All them are specified in such a way that smooth metric functions across the 
matching surface r = R are assured. For each set and any R, it is obtained 
M as a function of the scalar field amplitude. Taking as a reference the data 
set (a), for instance, it can be fitted the two others to it simply doing the 
appropriate rescaling in amplitude, as is showed figure [T] The limit mass for 
the background is R/2; this feature seems to be general. The limit mass to keep 
the incompressible static fluid is the well known Buchdhal mass limit AR/9. 
Therefore, if a black hole form it occurs by means of scalar radiation accretion. 
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Figure 1: Mass as a function of the normalized amplitude x — 'I'A^/ (6i?'^ +7rA^), 
for the exterior datum (a). If the datum is (b) the amplitude is A — ?> 13.46A; if 
the datum is (c) A — >■ ttA. For any choice of the initial datum and any choice of 
parameters r^, r;,, a and rg -which leave partially immersed the scalar field in 
the fluid distribution- the overlaping is true finding the appropriate rescaling 
in the amplitude. 



The exterior spacetime is determined integrating the field equations ( 29 1 



and (30). To evolve the scalar wave on the specified background we have the 



set of ordinary differential equations at r = for the scalar field ( 34 ) and ( 35 ) 



This system requires only initial conditions for $0j 'I'l and $2, which are fitted 
from the initial condition for the scalar field. Besides M and R are given as 
parameters, to define the matter distribution. 

The wave equation can be integrated following the null parallelogram method 



Note that, as is required by the used method, writing down (13) 

"ig.ur ~ [e^'^l - 2m/r)g,r],r = -[e'^(l - 2m/r)],,g/r, 



(68) 



where g — r$, we have to take care considering the RHS. To avoid the numerical 
derivative, that term is replaced by —At/, where 



o2/3 



A 



{47rr^(p - p) + 2rh/r] 



(69) 



A subtle issue of this model is the Killing propagator to get the energy con- 
servation as showed in figure [2j Because the spacetime is fixed as a background. 
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Figure 2: Energy conservation (multiplied by 10^) as a function of the Bondi 
time for the exterior /initial datum (a). The parameters and conditions of inte- 
gration are: M = 1.6921, = 5, A = 10 and 10^ radial points. The evolution 
corresponds to a relativistic case with F = 0.3232. The descending curve corre- 
sponds to A with M. calculated by Eq. ( 54 ) . The ascending curve corresponds 



to the energy radiated to infinity given by Eq. ( 58 ) with N calculated by means 
of Eq. (55). Thus, in accordance with Eq. (57), the horizontal curve represents 



the global conservation of energy. 



the right propagator is 

= e-^ff^p^. (70) 

Figure [S] displays the evolution of the initial data (a) ; it is evident that the 
Newman-Penrose constant is conserved. 



4.2 Scattering off a polytrope 

In the case of a static polytrope 



(71) 



where K is the polytropic constant and F is the adiabatic exponent, related 
with the adiabatic index n by F = 1 + 1/n, the pressure and the energy density 
vanish at r = i?. Again, we integrate numerically the system (10), (11) and 



(21) for any choice of K and F, tipically 100 and 2 respectively. However, it is 
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Figure 3: Decay of the scalar field (multiplied by 10"^) preserving the Newman- 
Penrose constant, that is, g^x\j+ = constant, any time of the evolution. The 
compactified coordinate x is related with r by means of r = 15i?a;/[8(l — x^)]. 
The initial datum, the parameters and conditions are the same as for figure [2] 



interesting to note that (21 1 written as 



P,r 

P + P 

can be integrated to get 



ln(l - 2m/r) + 2/3 



=2/3 



^^1/2(1 -2m/r)~i/2 



(1+p/p) 



n+l 



(72) 



(73) 



This result let us to conclude that no limit appears for the total mass distribution 
except the black hole mass limit R/2. It clearly connects the surface with the 
center of the distribution by means of e'^^°{l + Po/po)^'^^ = F^^^- 

In this model, a polytrope with a vacuum Schwarzschild exterior as a back- 
ground, the scalar field is scattered off and radiated to infinity. Figure [3] shows 
the energy conservation in this case. It should be stressed here that the scalar 
field gradient of the initial datum has to be zero at r — R and the scalar field 
itself partially immersed in the distribution. 
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Figure 4: Energy conservation (multiplied by 10'^) as a function of the Bondi 
time for the initial datum (c). In this case the matter fluid corresponds to a 
polytrope with po = 4 x 10^"^, K — 100 and F = 2. The parameters and condi- 
tions of integration are: M = 0.774, R = 11.512, A = 10"^, a ^ R, tq ^ R and 
10^ radial points; the evolution corresponds to a F = 0.865. The descending 
curve corresponds to A with A4 calculated by Eq. (54 1. The ascending curve 



corresponds to the energy radiated to infinity given by Eq. ( 58 ) with Af calcu- 
lated by means of Eq. (55). Thus, in accordance with Eq. (57 1, the horizontal 



curve represents the global conservation of energy. 



4.3 Quasinormal mode and late time tail decay 

For a Schwarzschild-Schwarzschild background an expected feature is showed 
when an initial compact support scalar field (b) is evolved. In this case the distri- 
bution of matter corresponds to an incompressible fluid, that is, a Schwarzschild 
interior spacetime; the exterior is a Schwarzschild vacuum. Figure [5] display the 
quasinormal mode ringing and the final tail decay. The energy conservation is 
displayed in figure [6j 



5 Concluding remarks 

We have used corner stones in the characteristic formulation of general and 
numerical relativity to present in this paper a framework which couples matter 
with radiation [37], [12]; [13] • The old point of view of Bondian observers are 
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Figure 5: Quasinormal mode and tail decaying as a function of the Bondi time 
for the initial datum (b). In this case the matter distribution corresponds 
to an incompressible fluid, that is, a Schwarzschild interior; the exterior is 
a Schwarzschild vacuum. The parameters and conditions of integration are: 
M = 1, i? = 2.857, A = 1, = 0, rf, = 10, and 10'^ radial points; the evolution 
corresponds to a F = 0.3. The late time behavior of the signal decays as an 
inverse power law. In the window of time [200, 245] the power of the tail decay 
is approximately —2.1, which corresponds to a 5% deviation from the expected 
value of-2 [35], [TD]. 



mistakenly considered as Eulerian in numerical relativity for the mathematical 
treatment of matter [?!] . 

Although we only explored toy models, we believe they deserve future at- 
tention at least as initial conditions. The partially immersed scalar field in the 
fluid is hidden. This mixture constitutes a fixed and well behaved background 
to evolve a cloned scalar field. The face value of model A is its usefulness as a 
playground to study asymptotic and matching regions. It is striking how the 
matching leads us to a connection between the mass distribution and the ampli- 
tude of the initial scalar field, and how the limit mass of R/2 does not depend 
apparently upon the initial data. The energy conservation and the Newman- 
Penrose constant tests give confidence on the possible physical consequences of 
this simple model. The other two models represent simple tests which are in 
agreement with expectations. 

We stress that this paper does not develop a general code, it just reports how 
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Figure 6: Energy conservation as a function of the Bondi time for the initial 
datum (b) . In this case the matter distribution corresponds to an incompressible 
fluid, that is, a Schwarzschild interior; the exterior is a Schwarzschild vacuum. 
The parameters and conditions of integration are: M — 1, R — 2.857, A = 1, 
Ta — 0, Tb — 10, and 10"^ radial points; the evolution corresponds to a F = 0.3. 
The descending curve corresponds to A with M calculated by Eq. (54). The 
ascending curve corresponds to the energy radiated to infinity given by Eq. ( 58 1 
with J\f calculated by means of Eq. (55). Thus, in accordance with Eq. (57), 
the horizontal curve represents the global conservation of energy. 



Bondian observers are plausible in Numerical Relativity. Numerics is simplified 
to show how these observers are consistent with previous results. The underlay- 
ing problem, gravitational radiation coupled to matter, is in fact difficult. We 
show that even in the spherical case we have to consider special regions, center 
and boundary surface carefully as well. We are currently developing a general 
code based on section 3, which considers regularization and matching regions in 
a clearer manner. 

The most simple axisymmetric case from the Bondian point of view is out of 
scope of the present investigation. There is a unique geometrical way to define 
Bondian observers in the absence of spherical symmetry: i) choose a tetrad 
to get the local Minkowskian frame; ii) make a double Lorentz boost to go to 
the frame comoving with the fluid. In this sense could be interesting to revise 
results from [2] to consider gravitational radiation considering explicitly Bondian 
frames and regularization on special regions such as the axis of symmetry and 
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boundary surface. 

An advanced study of incompressible-like fluids and polytropic matter indi- 
cates the relevance of central equation of state at the center of the star. This 
motivated us to prepare a hydrodynamic solver which will be reported elsewhere 
for an adiabatic situation. A more general study of matter coupled to radiation 
within this framework is in progress. 
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